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Abstract 

Isoradial dimcr models were introduced in [12] - they consist of dimer models 
whose underlying graph satisfies a simple geometric condition, and whose weight 
function is chosen accordingly. In this paper, we prove a conjecture of [12], namely 
that for periodic isoradial dimer models, the growth rate of the toroidal partition 
function has a simple explicit formula involving the local geometry of the graph 
only. This is a surprising feature of periodic isoradial dimer models, which does 
not hold in the general periodic dimer case [14]. 

1 Introduction 

In this paper, we solve a conjecture of [12], namely that the growth rate of the partition 
function of periodic isoradial dimer models can be expressed explicitly, using only the 
local geometry of the underlying graph. This is a surprising feature of isoradial dimer 
models, which does not hold in the case of general periodic dimer models. Indeed, in 
the general periodic case, Kenyon, Okounkov and Sheffield [14] obtain an expression 
for the growth rate of the partition function involving elliptic integrals, thus making it 
hard to do explicit computations. An example of application of our result yields a two 
line proof of Kasteleyn's celebrated result for the growth rate of the dimer partition 
function of the quadratic lattice Z 2 [7, 8], see Section 1.5. In order to state our result, 
let us precisely describe the setting. 

1.1 Dimer model 

The dimer model belongs to the field of statistical mechanics, and represents the 
adsorption of diatomic molecules on the surface of a crystal, it is defined in the following 
way. The surface of the crystal is modeled by a graph G. We assume that G is simple, 
infinite, simply connected (i.e. it is the one skeleton of a simply connected union of 
faces) and that its vertices are of degree > 3. A dimer configuration of G is a 
perfect matching of G, that is a subset of edges M, such that every vertex of G 
touches a unique edge of M\ refer to Figure 1 for an example in the case where G is a 
finite subgraph of Z 2 . Let us denote by M.{G) the set of perfect matchings of G. 
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Figure 1: A dimer configuration of a finite subgraph of Z 2 . 



Consider a simply connected finite subgraph Go of G. Then, dimer configurations of 
Go are chosen with respect to the Boltzmann measure, defined in the following way. 
Assume that a positive weight function v is assigned to edges of G, that is, each edge 
e of G has a weight v(e). Then, each dimer configuration M of Go has an energy 
£{M) = — X^eeM 1°6 v ( e )- The probability of occurrence of the dimer configuration M 
chosen with respect to the Boltzmann measure no, is given by: 

where Z(Gq,v) = Yl<MeM(G ) lleeM K e ) * s * ne normalizing constant, known as the 
partition function. 

The first step in the study of the dimer model is the computation of the partition 
function. Indeed, the latter yields precious information about the global behavior of 
the system [1, 2, 14], whose understanding is the goal of statistical mechanics. Kasteleyn 
[7, 8], and independently Temperley & Fisher [20], laid the ground stone for the study 
of the partition function. Consider Go as above. Then they give an explicit formula 
for Z(Go,^) as the square root of the determinant of a matrix Ko, also known as 
a Kasteleyn matrix. Loosely stated, Ko is the weighted adjacency matrix of Go 
(weighted by the function u), where minus signs are added to coefficients in a suitable 
way, see Section 2.2 for details. 

Since edges of dimer configurations represent diatomic molecules, the next step is to 
understand the partition function of infinite graphs. The relevant question then, is the 
computation of its growth rate. More precisely, if {G n } is an exhaustion of G by finite 



graphs, the goal is to compute: 

C= lim \\ogZ{G n ,v). (1) 

n— >oo n 

Note that this limit is not universal, in the sense that it strongly depends on the choice 
of exhaustion {G n }. Moreover, in the case of exhaustions by planar graphs, this limit 
is hard to compute because of the lack of symmetry which reflects in the Kasteleyn 
matrix. This is the reason why computations are done on toroidal graphs, as explained 
in the next section. 

1.2 Toroidal dimer models 

Consider an infinite graph G as above, and suppose that G is bipartite, i.e that it admits 
a bipartite coloring of its vertices. Assume moreover that there exists a bi-dimensional 
lattice A, such that G is doubly A-periodic, that is the graph G as well as its vertex 
coloring are bipartite. Define G n to be the toroidal graph G/An, then {G n } is a natural 
exhaustion of the graph G. 

In the case where G is the square lattice Z 2 , Kasteleyn [7, 8] gives an explicit expression 
for Z(G n ,u) as a linear combination of the determinants of four Kasteleyn matrices. 
Tesler [21] generalizes this result to any graph G satisfying the above, see Section 2.2 
for details (he actually generalizes it to graphs embedded on genus g surfaces). In [14], 
Kenyon, Okounkov and Sheffield give an explicit expression for the limit C of equation 
(1) in the case of the natural toroidal exhaustion {G n } of G. For general periodic, 
bipartite graphs G, this limit depends on the combinatorics of the model and involves 
elliptic integrals. 

In the next section, we define isoradial dimer models - the sub-family of dimer models 
for which the expression for the growth rate C of (1) becomes surprisingly simple. 

1.3 Isoradial dimer models 

Isoradial dimer models were introduced in [12]. They are dimer models on graphs G 
satisfying a geometric condition called isoradiality. This notion first appeared in [5], 
see also [16], and is defined as follows. All faces of G are inscribable in a circle, and 
all circumcircles have the same radius, moreover, all circumcenters of the faces are 
contained in the closure of the faces. The common radius is taken to be 1. An isoradial 
embedding of the dual graph G* is obtained by sending dual vertices to the center of 
the corresponding faces. 

Recall that the energy of a dimer configuration depends on the weights assigned to 
edges of G. In the case of isoradial graphs, one considers a specific weight function 
v, called the critical weight function, defined as follows [12]. To each edge e of G 
corresponds a unit side-length rhombus -R(e) whose vertices are the vertices of e and of 
its dual edge (R(e) may be degenerate). Let R = U eg G-R(e). Then, define v(e) = 2sin#, 
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where 29 is the angle of the rhombus R(e) at the vertex it has in common with e; 9 is 
called the rhombus angle of the edge e. Note that v(e) is the length of e* the dual 
edge of e. 

When the graph considered for a dimer model is isoradial, and when the energy of 
configurations is determined by the critical weight function, we speak of an isoradial 
dimer model. 

1.4 Result 

Consider an isoradial dimer model. Assume moreover that G is bipartite and doubly 
A-periodic. As above consider the natural exhaustion {G n } of G by toroidal graphs. 
Then, Theorem 1 below gives an explicit expression for the growth rate of the partition 
function of the exhaustion {G n } as a function of the local geometry of the graph. 

Theorem 1 

1 m (9 1 \ 

hm — log Z(G n , v) = V U log 2 sin 9> + -L(9i) , (2) 
n^oo n z r— f V vr it J 

where, 0\,...,6 m are the rhombus angles of the edges e±,...,e m of G\, and L is 
Lobachevsky 's function, L(x) = — log 2 sintdt. 

• The right hand side of (2) is (up to a factor 2), the expression Kenyon obtained 
for what he calls the log of the normalized determinant of the Dirac operator [12]. 
He conjectured it to be the right limit for the growth rate of the partition function 
of an appropriate exhaustion of G. Theorem 1 states that this is the right limit 
in the case of the natural exhaustion {G n } of G by toroidal graphs. 

• The surprising feature of Theorem 1 is that the growth rate of the partition 
function is expressed using geometric information of the fundamental domain G\ 
only, no combinatorics is involved. This is contrast with the result of [14] for 
general bipartite periodic graphs. Note that, although the result of [14] is true in 
a more general setting, there seems to be no simple way of re-deriving Theorem 
1 from it. Let us also mention that we use the result of [14] in order to have a 
priori existence of the limit C of equation (1). 

1.5 Example: Kasteleyn's computation 

In [7, 8], Kasteleyn skillfully computes the toroidal partition function for G = 1? . As an 
example of application of Theorem 1, let us rederive his result in a two line calculation. 
Note that the following computation was already mentioned in [12], as a support for 
the conjecture proved in Theorem 1. When G = 1? , the critical weight function is 
v = y/2, see Figure 2 (left). Moreover, consider the fundamental domain Gi = G/A as 
in Figure 2 (right). 
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8 = 71/4 
v (9) = 4l 



Figure 2: Critical weight function for I? (left). Fundamental domain G\ (right). 
Then, a direct computation using Theorem 1 yields: 



Acknowledgments: we would like to thank Richard Kenyon for proposing the question 
solved in this paper, and for fruitful discussions. 

2 Proof of Theorem 1 

In the whole of this section, we let G be an infinite, bipartite, A-periodic, isoradial 
graph. We assume that edges of G are assigned the critical weight function. Moreover, 
W denotes the set of white vertices of G, B the set of black ones, and G n the toroidal 
graph G/An. 

Before giving the proof of Theorem 1, let us precisely describe the conjecture of [12]. 
In [12], Kenyon introduces the Dirac operator K indexed by vertices of G (we refer the 
reader to [12] for the definition of K, and to Section 2.1 for the definition of the real 
Dirac operator K which is closely related). He also defines its inverse if -1 , for which 
he proves existence and uniqueness. One of his beautiful result is an explicit expression 
for the log of the normalized determinant "log deti K" of the Dirac operator K, defined 
by the following PDE: 




where K, is Catalan's constant. Note that the |log2 factor differs from Kasteleyn's 
computation, because we have weights y/2 on the edges instead of weights 1. 



d log deti K 



1 



K-\b,w), 



(3) 



dK(w,b) 



\V{Gy)\ 



where w G W, b £ B are adjacent vertices, |V(Gi)| is the number of vertices of G\. 
Adding some initial condition allows him to solve the PDE, and obtain the explicit 
expression for logdeti-RT, which is the right hand side of equation (2). How is this 
related to the dimer model? The Dirac operator K is closely related to an infinite 
Kasteleyn matrix - to have an actual infinite Kastcleyn matrix, one needs to work with 
the real Dirac operator K (see Section 2.1 below), obtained from the Dirac operator by 
a gauge transformation. If Go is a finite subgraph of G, and Ko the restriction of K to 
the vertices of Go, then logdeti Ko can be defined naturally as \y(Q j\ logdet Ko, which 
is the normalized log of the dimer partition function of Go- Hence, the conjecture is to 
interpret (3) as the limit of the dimer partition function on some appropriate exhaustion 
of G. The problem lies in proving the existence of such limits, leading us to work with 
toroidal exhaustions. Note that part of the proof is close to [12], since it consists in 
solving the PDE (3) for operators restricted to subgraphs of the exhaustion. 

The structure of Section 2 is as follows. Section 2.1 consists in the definition of the real 
Dirac operator K and its inverse K _1 . The operator K was introduced in [3], see also 
[13]. Note that features of K are closely related to those of the Dirac operator K of 
[12]. In Section 2.2, we state the explicit expression of [8, 21] for the toroidal partition 
function. Section 2.3 uses Sections 2.1 and 2.2, and consists in the proof of Theorem 1. 

2.1 Real Dirac operator 

The real Dirac operator K is obtained from the Dirac operator K of [12] by a gauge 
transformation. Both operators are represented by weighted infinite adjacency matrices 
indexed by vertices of G. For K the edges are unoriented and weighted by the critical 
weight function times a complex number of modulus 1. For K, edges are oriented 
with a clockwise odd orientation (see below), and are weighted by the critical weight 
function. Both weight functions yield the same Boltzmann measure on finite simply 
connected sub-graphs of G, but, and this is the reason why the real Dirac operator K 
is introduced, these weights do not yield the same probability distribution on toroidal 
subgraphs. Note that the real Dirac operator was already used in [3], see also [13]. 

2.1.1 Definition 

Following Kasteleyn [8], let us define clockwise odd orientations on edges of G. An 
elementary cycle C of G is said to be clockwise odd if, when traveling cw (clockwise) 
around the edges of C, the number of co-oriented edges is odd. Note that since G is 
bipartite, the number of contra-oriented edges is also odd. Then, an orientation of the 
edges of G is defined to be clockwise odd if all elementary cycles are clockwise odd. 
In [8], Kasteleyn shows that, for planar simply connected graphs, such an orientation 
always exists. 

Consider a clockwise odd orientation of the edges of G. Define K to be the infinite 
adjacency matrix of the graph G, weighted by the critical weight function v. That is, 
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if v\ and V2 are not adjacent, K(v±, V2) = 0. If w G W and b £ B are adjacent vertices, 
then K(w,b) = —K(b,w) = (—l) l ( w > b ^v{wb), where I( w ,b) = if the edge wb is oriented 
from w to b, and 1 if it is oriented from b to w. The infinite matrix K defines the real 
Dirac operator K: C v ^ -»• C y ( G ), by 

(k/)(«) = x;kk «)/(«). 

The matrix K is also called a Kasteleyn matrix for the underlying dimer model. 
2.1.2 Inverse real Dirac operator 

Similarly to the definition of the inverse Dirac operator K~ l [12], the inverse real 
Dirac operator K _1 is defined to be the unique operator satisfying: 

1. KK- 1 = Id, 

2. w) -► 0, when \b — w I — > 00. 

The rational functions f wx (z) are the analog of the rational functions f wv (z) of [12], and 
are defined in the following way. Let w 6 W, and let x £ B (resp. x £ W)\ consider 
the edge-path w = w±, 61, . . . , w^, b^ = x (resp. w = wi, b\, . . . , w^, bk,Wk+i = x) of 
G* from w to x. Let R(wjbj) be the rhombus associated to the edge Wjbj, and denote 
by Wj,Xj,bj,yj its vertices in cclw (counterclockwise) order; e tOL i is the complex vector 
t/j—Wj, and e i ^ j is the complex vector Xj—Wj. In a similar way, denote by Wj+i , x'j , bj,y'j 

the vertices of the rhombus R(wj+\bj) in cclw order, then e ta j is the complex vector 
y'j — Wj + i, and e^i is the complex vector x'j — Wj + ±. The function f wx (z) is defined 
inductively along the path, 

fww{z) = 1) 

IT . a -j+P-j 

(-l^K^'V^V^ 
fw6, (z) = ftuTO,- (2) — Zi^TTT. 



(z - e ia i){z - eW) 
f WWj+1 (z) = f^Wt-lj^^Ve-V^-e'^tz-e^). 



By Lemma 9 of [3], the function f wx is well defined. Moreover, Lemma 10 of [3] gives 
an explicit expression for K _1 (it is the analog of Theorem 4.2 of [12]): 

Lemma 2 [3] The inverse real Dirac operator is given by 

K" 1 (b, w) = J f w b(z) log z dz, (4) 

where C is a closed contour surrounding cclw the part of the circle {e td \ 9 € [0q — vr + 
A, #0 + 7T — A]}, which contains all the poles off w b, and with the origin in its exterior. 



7 



2.2 Partition function of the torus 

In this section, we give the explicit formula of Kasteleyn [8] (in the Z 2 case), and Tesler 
[21] (in the general periodic case) for the toroidal partition function Z(G n ,u). Note 
that this result was already described in [3], we nevertheless choose to repeat it here 
since its understanding is important for the proof of Theorem 1. 

Let us first orient the edges of G. Consider the graph G\, it is a bipartite graph on the 
torus. Fix a reference matching Mq of G\. For every other perfect matching M of G±, 
consider the superposition M U M of M and Mq, then M U Mq consists of doubled 
edges and cycles. Let us define four parity classes for perfect matchings M of G\\ (e,e) 
consists of perfect matchings M, for which cycles of M U Mq circle the torus an even 
number of times horizontally and vertically; (e,o) consists of perfect matchings M, for 
which cycles of M U Mq circle the torus and even number of times horizontally, and an 
odd number of times vertically; (o,e) and (o,o) are defined in a similar way. By Tesler 
[21], one can construct an orientation of the edges of Gi, so that the corresponding 
adjacency matrix K} has the following property: perfect matchings which belong to the 
same parity class have the same sign in the expansion of the determinant; moreover 
of the four parity classes, three have the same sign and one the opposite sign. By an 
appropriate choice of sign, we can make the (e,e) class have the plus sign in detK}, 
and the other three have minus sign. Consider a horizontal and a vertical cycle of 
the dual graph G* of G±. Then define K\ (resp. Kg) to be the matrix K{ where the 
sign of the coefficients corresponding to edges crossing the horizontal (resp. vertical) 
cycle is reversed; and define K4 to be the matrix K} where the sign of the coefficients 
corresponding to the edges crossing both cycles are reversed. By Kasteleyn [8] (in the 
I? case), and Tesler [21] (in the general case), we have the following, 

Z(G u v) = -(- det K\ + det + det + det K\). 

The orientation of the edges of G± defines a periodic orientation of the graph G. For 
every n, consider the graph G n , and the four matrices K™, KJ? , K3 , K4 defined as above. 
These matrices are called the Kasteleyn matrices of the graph G n . 

Theorem 3 [8, 21] 

Z(G n , v) = \(- det K™ + det K% + det K™ + det KJ). 

The orientation defined on the edges of the graph G is a clockwise odd orientation. Let K 
be the real Dirac operator indexed by the vertices of G, corresponding to this clockwise 
odd orientation. Note that except for edges crossing the horizontal and the vertical 
cycle, the coefficients of the Kasteleyn matrices K" and of the real Dirac operator agree 
on edges they have in common. 



8 



2.3 Proof 



Consider an orientation of the edges of G denned as in Section 2.2, and let K be the 
real Dirac operator indexed by vertices of G, corresponding to this orientation. Let 
K™, K2 , Kg , K2 be the Kasteleyn matrices of the graph G n . The following computations 
are inspired form [12]. By linear algebra, for every I = 1, . . . , 4, and for every edge W{bi 
of G n , we have 

9d6tK " = (det ^((K?)- 1 ^, «,,))■ 



d(K?(wi,bi)) 

Denote by e± = w±bi, . . . ,e m = w m b m the edges of G±, and let 9i,...,6 m be the 
corresponding rhombus angles. Since the graph G n is invariant under A-translates, 
we know that for every A-translate w\b\ of the edge lOj&j, the coefficient K™(w*,6*) = 
±Kg(wi, bi). The minus sign only occurs when I = 2, 3, 4, (recall that in the definition 
of K?> , . . . , K", the sign of the entries which cross the horizontal and/or the vertical 
cycle of G* is reversed), but then (K™) _1 (6*, w\) = ±(K™) _1 (6j, Wi). Hence, for every 
t = 1, . . . , 4, and for every i = 1, . . . , m, 

ddetK™ ^ ddetK™ d(K n t (w\,b^) 

w\o\ translates or Wibi 

= n^WWrH^VH)) ^^ , (5) 

where the edg 6S 61 , . . . , G m do not cross the horizontal and the vertical path of the dual 
graph G* of G n . Define the function ip n . . . , 9 m ) by 

<Pn(0l,---,0m) = \\ogZ(G n ,v). 

n z 

By Theorem 3, we have Z(G n , v) = \{- det KJ + det K£ + det + det K£). Moreover, 
for every i = 1,...,4, K"(itfj,&j) = K(wi,bi), so that using (5), we obtain for every 
i = 1, . . . , m, 



(6»i,...,6> m ) - ^ — J ^——-----(K!) (b I)TOl ) + ^ 2Z(Gri i/) (K,) (&,,«>,) 



(6) 



The next part of the argument can be found in [9]. The second bracket of equation (6) is 
a weighted average of the four quantities (K") _1 (6j, Wi), with weights ±| det Kg/Z(G n , v). 
These weights are all in the interval (—1/2, 1/2) since for every I = 1, . . . , 4, 2Z(G n , v) > 
|detK™|. Indeed, Z(G n ,u) counts the weighted sum of dimer configurations of G n , 
whereas | det K"| counts some configurations with negative sign. Since the weights sum 
to 1, the weighted average converges to the same value as each (K™) _1 (6j, Wi). 

By Proposition 2 of [3], for every t = 1, . . . , 4, (K™) _1 (6j, Wi) converges to K~ 1 (6j, Wi) 
on a subsequence (rij) of n's. Hence, for every e > 0, there exists uq such that for 
n > uq, n £ (rij), equation (6) can be written as 

-^-(6>i,...,0 m ) = 7^ K (bi,Wi)+e. (7) 
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Let us compute the right hand side of (7), using the notations of Section 2.1.2. By- 
definition, K(wi,bi) = (-l/K.M 2 sin 9 i} hence 

= (-1)^)2008^. (8) 

OVi 

Moreover, by definition f Wi ^(z) = (—1) (^-W) ^ z _^ eiai ^ z _ el(3i ^ - So that using Lemma 2, 
and the Residue Theorem yields: 

K-Hk,m) = ' (-l)'.-,..)e^ / , fi' to, 

47rn J c (z — e ia *)(z — e llJ *) 



V e^ ai — e*^* gift— a 

_L (-IjW*) -A- (9) 
2vr v ; sin ft- V ; 



Combining equations (8) and (9), implies: 

dip. 



(01, . . . , m ) = — cotanft + e. 



00, - 



By [12], there is a continuous way to deform the graph G so that all rhombus angles 
tend to or ir/2. The same transformation can be applied to G n , for every n. Denote 
by Of the angle 0{ after such a deformation. Let M be the number of angles 6®, . . . , 9 m 
which are equal to tt/2. By the argument of [12], for n > no, n € (rij), we obtain 

m / 6 ■ 1 \ M 

<p n {0i,...,e m ) = Y, (-^28^ + -^)) - — i og 2 + ^„(e?,... + (io) 



Let us compute <p n {Q\i ■ ■ ■ , 9 m )- Suppose the above deformation is applied to the graph 
G, then edges corresponding to rhombus angles have weight 0, and those correspond- 
ing to rhombus angles ir/2 have weight 2. Removing the weight edges, the deformed 
graph consists of independent copies of Z. Applying the deformation to G±, we obtain 
graphs Z/miZ, . . . , Z/m p Z, for even mi, . . . , m p , where each TLjmflL consists of weight 
2 edges. We can compute, 

nm-i nrrir) o M 

Z(G n , v) = (2^2) n . . . (2^2)™ = 2" ~2 np . 
Hence, for every e > 0, there exists n 1; such that for every n > we have 

^,...,0 = y!og2 + e. (11) 

Combining equations (10) and (11) yields, for every n > max{no,rai}, n G (n^), 

m (B 1 \ 

<^(0i, . . . , e m ) = V p log 2 sin 0, + -L{6i) + 2e. 

i=l V 7 
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This implies, 

m (Q 1 \ 

lim' n ^ ooV3ri (0 1 , . . . ,9 m ) = V - log 2 sin ^ + -L(^) . 

z — ' \ 7T 7T / 

1=1 v ' 

where the limit is taken on the subsequence (jij). Using Theorem 3.5 of [14], we deduce 
that <p n (0i,..., 9 m ) converges (as n — > oo) to the same limit as the above subsequence. 

□ 
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